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Abstract 



The ^-meson production cross section is measured for the first time at a sub- 
threshold energy of 1.93 AGeV in ^^Ni+^^Ni central collisions. The (/) data were 
obtained within the acceptance of the CDC/Barrcl subsystem of FOPI. For a sam- 
ple of 4.7 • 10^ central events, after background subtraction, 23 candidates were 
observed. Extensive GEANT simulations of the detector performance are shown 
in a thorough comparison to the real response, aiming at a good understanding 
of the apparatus and at a trustable determination of the efficiencies, production 
probability and possible systematic errors. A filter procedure is elaborated, which is 
meant to facilitate the comparison of any theoretical calculation or new data with 
the current ones. How to extrapolate the present value to a ^- meson cross section 
in 4:71 is also discussed. This result on pseudo-vector mesons can now be compared 
to existing experimental knowledge for the same reaction at the same incident en- 
ergy for various outgoing channels, K'^ and K~ included. A significant fraction (at 
least 20%) of the i^~-mesons is originating in the decay of the (f), supporting the 
statement that the two channels are strongly correlated. 

Key words: heavy ion collision, SIS, sub-threshold meson production, 
pseudo- vector mesons, (/)-meson 
PACS: 25.75.-q, 25.75.Dw 



1 Introduction 



Relativistic heavy ion collisions offer a unique possibility to probe experimen- 
tally hot and dense nuclear matter: at bombarding energies of 1 — 2 AGeV 
baryonic densities of 2 — 3 times the saturation value po at zero temperature 
are reached, however the condition lasts merely about 15 fm/c and the degree 
of equilibration established in such a short interval is open to questions [1]. 
A long standing goal of the study of central reactions in the beam energy 
range 1 — 10 AGeV is the determination of the compressibility k of nuclear 
matter at non-zero temperature and high density. Unfortunately, almost all 
the proposed observables sensitive to k have at some point been discovered 
unreliable. Many forms of collective motion (flow [2]) were thought to mani- 
fest simply the effect of the compressibility of nuclear matter, but were later 
proven to be strongly affected also by the momentum dependent part of the 
mean field so that the former one can not be simply disentangled [3,4,5] (in 
ref. [5] a reasonable estimate of 165 < k, < 220 MeV was derived from the 
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data). Recently elliptic flow has been proposed as a possible test quantity [6]. 
Pion production was believed [7,8,9] to be another good observable to extract 
k; but the high absorption in nuclear matter makes it susceptible not only 
to the compression phase and maximum density but also to the details of 
the expansion [10,11,12]. Finally, Aichelin and Ko [13] proposed sub-threshold 
K^-production as a promising probe: in fact strangeness is conserved in strong 
interactions and absorption is negligible in nuclear matter (strangeness ex- 
change is present but at equilibrium should be balanced) and hence the yield 
is insensitive to the expansion phase. Because of the Pauli blocking the di- 
rect nucleon-nucleon reaction, that could happen with the help of the Fermi 
momentum, is almost suppressed and at sub-threshold energies the dominant 
creation mechanism is a two step one. The conclusion has been confirmed ex- 
perimentally by the KaoS collaboration measuring the variation in the yield 
with the number of participants to the central coUision [14,15,16]. The pro- 
duction probability is then tightly bound to the maximum density reached in 
the compression phase. 

However such apparently clear situation is obscured by the possibility of in- 
medium effects on the K~^, K~ and in general meson properties. Despite the 
small alteration of the K'^ mass in the nuclear environment, the experimen- 
tal determination of the EOS is rendered difficult [17]. Studying the relation 
between system size and i^'"'"-production with high quality data, in order to 
distinguish between the in-medium and fire-ball density (i.e. EOS) contribu- 
tions, has accumulated an evidence in favor of a soft equation of state of 
nuclear matter {k ^ 200 MeV) [18]. 

Nevertheless the subject of the possible infiuences of the surrounding nuclear 
medium on the meson masses is interesting by its own, with important impli- 
cations in astrophysics [19,20,21,22]. For example, eventually the conden- 
sation in neutron stars [23,24,25] could lead to a softening of the equation of 
state of nuclear matter above 3 times the saturation density Pq. The evolution 
of supernova explosions is than altered causing a core with 1.5 — 2 solar masses 
to collapse into a black hole rather than to form a neutron star [20,21,22]. 

In particular it is in general well accepted that in nuclear matter at equi- 
librium increasing temperature and/or density will lower the so called quark 
condensate {qq) thus partially restoring the chiral symmetry. However, while 
the prediction of the consequences on pseudo-scalar mesons (like the and 
K') is robust to different model frameworks [26], the situation appears dis- 
puted for pseudo- vector mesons (like the (p) more exposed to problems in the 
theoretical description [27,28,29,30,31,32]. Additionally the 0-meson yield is 
coupled to the sum of the and K~ mass. If, as it is currently hypothe- 
sized, the reduction of the latter can be as high as ~ 100 MeV for p = 3 po, 
the corresponding decay width in the nuclear medium is diminished increasing 
the absorption probability. Actually Pisarski and independently Shuryak and 
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collaborators [33,34,35] originally proposed to study the kaon in-mcdium prop- 
erties from an enhancement or suppression of the 0-production probability in 
heavy ion collisions, implicitly considering the modification of the (p mass itself 
less important. Hartnack, Oeschler and Aichehn have suggested that the K~ 
yield is almost insensitive to its own in-medium mass, but is instead strongly 
linked to the one by the secondary interaction K'B ^ A(S)7r [36]. All 
these arguments indicate the strong interconnection between the , and 
0-channels and assess the impossibility of studying each of them separately. 
The present result is the missing piece to complete the already available knowl- 
edge on and X^-mesons [37,38,39] for the same reaction and beam-energy. 

Besides, due to a; — mixing, the 0-meson wave function is dominated by 
the ss component and so its propagation in nuclear matter is interesting for 
complementing the informations gained from hadrons with open strangeness 
{K+,K-,K%K) [40,37,41,39,42]. 

Recently the parameterization of the NN in-medium cross section has been 
shown to be the major element of uncertainty for the calculation of the charged 
kaons production probability. As a matter of fact the disagreement between 
the Giessen group [43] and Nantes group [4,44] transport codes, quite dis- 
similar on a physical basis, was found to be marginal, under such respect, 
when a common parameterization is supplied [36] . It is also comparable to the 
predicted variation between the in-medium and no in-medium effect scenarios. 

The long awaited measurement of both the elementary processes pp —> K^X 
[45,46] and pp — > (ppp [47,48] near threshold allows to nail down in a crit- 
ical region the freedom in the cross section parameterizations necessary for 
theoretical models, consistently reducing ambiguities. 

In this work the data from FOPI restricted to the CDC/Barrel subsystem on 
0-meson yield in ^^Ni-|-^^Ni at 1.93 AGeV will be introduced. Prom within the 
limited experimental acceptance, a 47r extrapolation can exclusively be done 
once the shape and characteristics of the source are assumed. The channels 
and thresholds for (/)-meson creation are not too different from that of a K~ 
pair and, as recalled, they are strongly related: our best estimate of the (j)/ K~- 
ratio for central collisions will be discussed. A significant part of the K" , 
at least 20%, seems to arise from the decay of the 0. Preliminary hints of 
sub-threshold 0-meson production have already been reported using both the 
CDC/Barrel [49] and the Helitron/Forward Wall [50] combinations, here the 
emphasis is on the extraction of a quantitative estimate of their occurrence 
probability in central collision. 
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2 Data analysis 



2.1 The FOPI detector 



The FOPI detector is a modular system for fixed target experiments at the 
SIS accelerator of GSI-Darmstadt. The four subdetectors (CDC, Barrel, He- 
litron and Forward-Wall [51,52]) are coupled two by two (CDC/Barrel and 
Helitron/Forward-Wall) to provide momentum, energy-loss and time of flight 
information in the polar angles 39° < i?iab < 133° and 7° < i?iab < 28°, re- 
spectively. The coverage in the azimuthal angle (/? is almost complete, except 
for the Plastic Barrel. The reference polar coordinate system of the current 
work has the beam direction (i.e. z-axis) as opening direction and the target 
as origin for the polar angle 'j?iab- The azimuthal angle (p lies in the plane 
orthogonal to the beam axis containing the target position (i.e. xy-plane). 

The Barrel and the CDC are operated in a magnetic field of 0.6 T produced 
by a super-conducting coil. The Helitron is a vector chamber and together 
with the Forward- Wall is placed in the fringe field of the solenoid in the for- 
ward hemisphere. In the ongoing analysis the CDC/Barrel subsystem alone is 
involved and in what follows we focus on it. 

The CDC is a drift chamber of the Jet type [53,54,55,56,57,58] divided in 16 
sectors with projective geometry for the forward edge (inner radius 20.7 cm 
and length 75 cm, outer radius 80.1 cm and length 192.5 cm, see table 1) 
each equipped with 60 sense wires in the median plane. Such choice has the 
advantage of allowing a more uniform charge collection and multiplication field 
together with a reduction of the drift length by a factor of 2, still keeping an 
equal amount of electronic channels. However, the effective readout comes on 
cost of the left-right ambiguity for every hit that for tracking purposes needs 
to be considered in a pair of spatial positions of which just one corresponds to 
the correct localization. To partly remedy the problem the sectors are tilted 
by 8° with respect to the ideal radial symmetry around the target so that only 
one of the two mirror images points to the latter. In addition, the sense wires 
are staggered [53,54,55,56,57,58] by ± 200 fim about the nominal sense wire 
plane giving an independent criterion for the distinction of real and mirror 
tracks. The gas mixture chosen is 88% Ar, 2% CH4 and 10% iso-butane at 
atmospheric pressure. 

In Jet chambers the transverse position of the hit in the xy-plane is deduced 
from the cell number and measured drift time, while that along the ^-direction 
by the charge division between the ends of the wire: the physical basis is quite 
different for the two coordinates and also the physical mechanisms imposing 
the actual limits to the resolutions are not the same; as a consequence the first 
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is known much better than the second. Such feature is taken into account in 
the processing algorithm where the hits candidate to form a track are identified 
initially in the xy-projection of the event [54] requiring the distance from the 
nearest neighbor to be less than ^xyhit (see table 2). In the successive step 
the position along z of the hits pertaining to a track is examined and all those 
further apart from the closest one of i^zhit (see table 2) are rejected as not 
belonging to it. For each particle the transverse momentum component pt over 
charge Z (in units of the elementary charge) can be obtained by the inverse 
k of the radius of the arc of circle fitted through the xy-projection of the 
trajectory, as 

Pt _ 0.3 B GeV/c T'^m'^ _ 0.18 GeV/c m~^ 

Z ~ k ~ k ' ^' 



where in the last expression we substituted B = 0.6 Tesla for the magnetic 
field of FOPI. Prom the angle ■j?iab, dehvered by the straight line fitted to 
the ^-coordinate of the hits, the total momentum p over charge Z and its 
components are found 

^ = (2) 
Z sin i?i,b ^ ' 



The error on p is 



1 



(3) 



where the two autonomous contributions are apparent: a^, coming from the 
determination of the radius of curvature of the track in the transverse plane 
and, cTijj^j^, from the reconstruction of the emission angle ■j?iab- In the CDC 
the physical origins of and o^^^^ are quite dissimilar as we have already 
discussed. The condition 2/k > Rmm for the trajectory to enter the CDC can 
be reexpressed with eq. (1) as pt > 0.019 GeV/c. 

The particle identification is based on a double technique allowing a redun- 
dancy very powerful to reject the background. Once the momentum is known, 
the charge over mass ratio can be determined (and hence the species assigned) 
if the velocity is measured. 

Exploiting the dependence of the mean energy-loss on the velocity, and not on 

the momentum, of the particle sweeping through the detection medium (i.e. 
the CDC gas), using an empirical parameterization, a mass named Cmass 
can be extracted. In fig. 1 (upper two panels) the truncated mean energy- 
loss [54,56] (after correction for the ■& angle of incidence in each cell evaluated 
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from the track fit) is shown versus momentum: the branches corresponding to 
TT, K and p are discernible and separable. Since when the Barrel is reached 
the track length inside the CDC and hence the number of sampling points for 
the energy-loss distribution are constrained, only such cases are considered 
in fig. 1, besides it will be the final requirement of the present analysis. The 
continuous curves, also drawn, refer to the empirical parameterization defining 
the regions where the various species can be identified. 

An uncorrelated information on the mass is derived for the ones that survive up 
to the Barrel by the time of flight and the length of the extrapolated trajectory 
fitted to the CDC hits (a condition is imposed on the difference in <^-angle 
and in ^-coordinate between the measured and expected impact position 
in the Barrel, see table 2). In this case the velocity is determined directly. The 
Barrel consists of 180 bars of scintillator of width 3.25 cm, depth 4.05 cm and 
length 240 cm arranged in 30 modules (the space for 2 further ones is occupied 
by the CDC support structure) disposed in a cylindrical fashion at a radius 
of i?Bar = 111-5 cm coaxially with the z-axis and covering the -i^iab angles 
reported in table 1. Due to the missing modules and to the inter-module and 
inter-bar spacings also the azimuthal coverage is incomplete. The scintillators 
are readout with photo-multipliers located outside the FOPI magnet through 
long light guides. The condition 2/k > i?Bar for the Barrel to be reached can 
be reexpressed with eq. (1) aspt > 0.10 GeV/c. In table 1 we have summarized 
the acceptances of the CDC and Barrel subdetectors of FOPI. In fig. 1 (lower 
two panels) the velocity is displayed versus momentum and again the branches 
corresponding to tt, K and p are distinguishable. The continuous curves are 
the well known relativistic expression in the nominal cases. A mass named 
Bmass is calculated as 



The Cmass is limited at high momentum from the relativistic rise of the 
energy-loss as a function of p: above ^ 0.6 GeV/c the pion and kaon region 
merge and contamination is unavoidable (see fig. 1 a and b). 

Also Bmass is hmited, for kaons, below 0.6 GeV/c by the finite time res- 
olution of the Barrel (~ 400 ps) joined to the divergence in the relativistic 
momentum- velocity relation (see fig. 1 c). 

The advantage of the redundant simultaneous measurement of both Cmass 
and Bmass, is illustrated in fig. 2 representing one mass versus the other: 

a structure corresponding to the is clearly visible, the K~ are less well 
separated owing to their small abundance. To perceive an idea of the difficulty 
of such task, let us remind that in the acceptance for 1 proton there are roughly 
10~^ pions, fti 10~^ K'^ and 10~^ K~. The superior quality of Bmass 




(4) 
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over Cmass can be deduced from the deformation of the distribution along the 
latter. For this reason they will not be treated alike in the final identification 
cuts. 



2.2 Event selection 



The online event centrality selection is based on the detected charged particle 
multiphcity in the Forward- Wall, M, obtained with fast analog electronics [51]. 
A minimum bias trigger is defined when more than 4 are found, and a central 
one when M > Mc- During event acquisition, the quantity Mc is set so that 
the ratio between the minimum bias and central triggers corresponds to about 
10% [51]. In particular we focus on the 7.7-10® central collisions collected with 
the FOPI detector in 1995 for the ^^Ni+^^Ni reaction at 1.93 AGeV beam 
kinetic energy. After imposing constraints on the energy deposited in the start 
counter, for the purpose of optimizing the time resolution, and on the primary 
interaction vertex position in both the transverse xy-plane and ^-direction, 
for the purpose of reducing background from interactions not produced in 
the target, a sample of 4.7 • 10^ events is left. In the more accurate offline 
reprocessing a certain jitter is observed in the multiplicity threshold for the 
central trigger caused by the non-ideal behavior of the analog electronics with 
which it is implemented (see fig. 3). The similarity of a selection on minimum 
multiplicity M to an upper limit on impact parameter b{M) is straightforward 
in the sharp cut-off hypothesis assuming a monotonic relation between M and 
b. Adopting a simple geometrical model the central collision condition cited 
through this work can be seen to be equivalent to < 6 <~ 3.3 fm or to the 
most central (Tc ~ 350 mb (i.e. 12%) of the geometrical reaction cross section 
(Tgeo ~ 2.9 b [38]. A rather simple comparison is possible for the yield per 
central event between our experimental finding and transport codes once the 
the prediction is studied as a function of centrality. For the theoretical result 
-fth(^) an appropriate mean (Pth) over b in the range < 6 < 3.3 fm is to be 
computed taking as a weight the ratio of the differential reaction cross section 
dan/ db to the integrated one cxc 




max 








where 6, 



max 



is 3.3 fm in the case of FOPI. 
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2.3 ^ -meson identification 



The 0-meson is identified through its decay channel in a K'^ and a K~ 
{BR — 49.1%). The charged kaons are pre-selected mainly with a window- 
in both Cmass and Bmass and an upper momentum limit. A data sub- 
sample is generated with kinematic and other quantities. For calculating the 
invariant mass of the couple, the common nominal value of the two species 
is used. The combinatorial background is reconstructed via the event-mixing 
technique [59,60] : at the same time an independent data sub-sample is gener- 
ated regrouping in turn each kaon from the current event together with the 
opposite charged partner from a stack of depth 20 of previous ones. The latter 
are sorted with respect to the amount of particles detected in the Forward- 
Wall (i.e. impact parameter, see sec. 2.2). The stack can be filled with central 
collisions where a K'^ or a K~ is pre-identified or where both a and a K~ 
are pre-identified. The momentum vectors from the stack are rotated around 
the beam axis so that their reaction plane has the right current orientation. 
The population size in the event-mixing is easily made much larger (2 orders of 
magnitude) than in the data, in such fashion the statistical error contribution 
to the outcome for the ^-candidates is neghgible. 

To reject mirror or erroneous tracks, three general cuts to ensure good quality 
are imposed: on do (the minimum distance in the transverse xy-plane of the 
trajectory from the reconstructed event vertex position), on zq (the displace- 
ment from the target of the point where the helix crosses the 2;-axis) and on 
Uhit (number of hits): |(io| < 0.5 cm, |^o| < 10 cm and Uhus > 30 (see table 2). 
A difference between do and Zq should be noted: while the target thickness is 
small and the initial collision 2;-coordinate is well defined, the corresponding 
statement is not true for the beam spot transverse dimension and do refers to 
the xy-position extracted from the data. By the way, the acceptance windows 
on do and Zq also help to reduce the background of pions coming from sec- 
ondary vertices. All tracks fulfilling the do and Zo selections are then refitted 
forcing their origin in the event estimated first interaction point. This new 
helix is then the base for calculating the momentum and for performing the 
extrapolation to the Barrel necessary to evaluate the matching condition and 
the length of flight. In the present situation the procedure can not rise any 
problem since all the 0-mesons are decaying almost immediately outside the 
fireball (with cr = 44.5 fm and p = 1 GeV/c the mean length traveled is 43.7 
fm) and hence they are all well inside the target. It is of course not accom- 
plished in those studies, like on A-baryon production, in which detection of a 
vertex apart from the primary reaction one has to be exploited [40,42] . 

The final cuts (see table 2) piab < 0.6 GeV/c, Cmass > 0.37 GeV (for pion 
background suppression, see fig. 2 a and b), Cmass < 0.8 GeV (for proton 
background suppression, see fig. 2 a) and |Smass — 0.494| < 0.1 GeV are then 
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applied for particles of both charges to the two sub-samples (data and event- 
mixing) and the K~ invariant mass distributions are constructed (as shown 
in fig. 4). The normalization is provided integrating the spectra aside from the 
0-peak (minv > 1-06 GeV/c). The combinatorial background drawn with the 
lines enclosing one standard deviation in fig. 4 is obtained filling the mixing 
stack with events containing either a or a K~ ^ the other, represented with 
the standard error bars, when the mixing stack is filled from events containing 
a and a K~ . The first is much more regular and will be employed in the 
subsequent analysis. The second has still instabihties of non statistical origin 
indicating that, although the whole sample is just 30% smaller, the poorer 
class from which it is drawn (i.e. the same element of the stack is reused more 
often) makes the destruction of the correlations harder. It can be considered 
solely for estimating a kind of maximum systematic error due to self correlation 
induced by the limited population adopted for mixing (such source would not 
be important in an analysis with much more candidates). 

After subtraction, the number of events in the 0-peak (defined from the end- 
points at 1.01 and 1.03 GeV where the remainder of the signal and the back- 
ground becomes negative) with its statistical error is of 23 ± 7. In a previous 
analysis [61], done with older cahbrations, a result of 22 was found. Prom this 
difference and from the difference obtained subtracting both backgrounds a 
systematic error of ±2 or 9% is estimated. 



3 Simulation 

The (23 ± 7 ± 2) 0-mesons given in sec. 2.3, are of little utility for a com- 
parison with future experiments or with theoretical calculations. An accurate 
determination of the detection efficiency by simulation is needed to deduce a 
cross section within the phase space coverage (or in 47r). 

In the present section we will always assume that all (/)-mesons are decaying 
through the i^~-channel. The appropriate branching ratio will be included 
in sec. 4 for deducing the production probability. 

The overall efficiency etot to measure in our device such particle out of 47r can 
be splitted into two terms 

Ctot = Cmax Cdet ■, (6) 

where Cmax is defined as the efficiency of measuring a within the geometrical 
and momentum acceptance of a "perfect" apparatus plus the effect of kaon 
decay in flight during their path to the Barrel and edet contains the further 
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reduction introduced by detector limitations and by the constraints necessary 
for a clean identification (vertex, matching and mass cuts). 

Information on the (f) distribution in An can not be obtained in the current 
effort because the phase space region covered by the CDC/Barrel subsystem 
around mid-rapidity is small. For an extrapolation into the full solid angle of 
the cross section extracted in the acceptance, one has to rely upon a Att shape. 

The advantage of the factorization (6) is that while emax depends essentially on 
the initial phase space distribution in 4tt, Cdct does not, within a good approx- 
imation. Such case will be demonstrated to hold in a representative example 
in sec. 3.2 and 3.3 where the "effective" temperature of a one source thermal 
model is varied. Thus the comparison between our experiment and new ones or 
different calculations is greatly facilitated. For any elaborated study delivering 
the entire (f) population, it is sufficient to apply a simple "filter" procedure (see 
appendix A). This takes into account the detector geometry and momentum 
cut as well as the kaon decay, neglecting first all complicated details of the 
apparatus performance, which can be considered separately through the con- 
stant coefficient edet, determined once for all with the necessary knowledge of 
FOPI in sec. 3.3. 



3.1 The CDC/Barrel maximum acceptance 

In the study developed here the simple hypothesis of a thermal source con- 
trolled by only one parameter (i.e. the temperature T) is adopted. A more 
refined model would be the Siemens and Rasmussen [62] formula, contain- 
ing a temperature Tsr = 92 MeV and a collective expansion with velocity 
/5flow — 0.3. The whole idea behind this description of a central collision is 
that Tsr and /3flow a-re common for all hadrons and so their values can be 
supplied by a previous study done with the FOPI detector about p, d and n~ 
from the same ^^Ni+^^Ni reaction at the same beam-energy of 1.93 AGeV [63]. 
However, as long as just the 0-meson is concerned and within the necessary 
accuracies, the source can be approximated with a pure thermal one with 
T — 127 MeV. The influence of the dynamics can then be reabsorbed in the 
temperature T, that becomes distinct from the real one Tsr (and species as- 
sociated). Transport model calculations generally follow a similar trend: the 
shape of the phase space distributions can be considered thermal within the 
present statistical and systematic errors, but the effective temperatures are 
lower, typically 70 — 90 MeV. We can than in rough approximation reduce 
our uncertainty on the source to a thermal source where the temperature is 
unknown by a wide margin. 

To investigate the change in CDC and Barrel acceptances with temperature 
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T in the interval 130 - 70 McV. four GEANT simulations (see table 3) were 
undertaken. However the detector response description and event reconstruc- 
tion stages are not involved in the extraction of such estimates and solely the 
FOPI geometry implementation and particle transport part of the code were 
exploited. The results are summarized in table 3 with and without kaon decay 
enabled. The column corresponding to emax is the third from the end: a total 
excursion with temperature by a factor of ~ 4 is found. This clearly illus- 
trates the impossibility of giving one cumulative number and the necessity of 
the decomposition (6) in order to allow a comparison with the predictions. 

It is important to realize that the smallness of emax derives from the localization 
of the phase space covered by the CDC/Barrel around the target rapidities. In 
fig. 5 a the pt/m versus normalized rapidity y plot is shown for the and K~- 
mesons reaching the Barrel and originating in the decay of the emitted by a 
thermal source with 130 MeV temperature. As can be seen, the relevant cuts 
are the Barrel forward ^9iab angle limit of 39.2° and the maximum momentum 
allowance of 0.6 GeV/c (see table 1). In fig. 5 b the pt/m^y representation is 
used for the primary quantities of the 0-mesons whose and K~ decayed 
into the range of acceptance. Its maximum is still well away from the mid- 
rapidity region. The dot-dashed curves, also in panel b, correspond to the 
most probable momentum for a relativistic Maxwell-Boltzmann emission law 
centered at y = 0. The portion of the diagram accessible with the CDC/Barrel 
sub detector combination is outside this line indicating that exclusively the tail 
of the distribution is probed. This becomes even more pronounced if the source 
temperature is lowered, indicating why emax is reduced under such conditions. 



3.2 The detection efficiency estimation from GEANT 

As already mentioned e^et contains the adjoined contribution to the efficiency, 
in the phase space acceptance and after kaon decay, of the detector behavior. 
For its determination a precise reproduction of the apparatus response is nec- 
essary. The current section is devoted to the discussion of the quality of the 
environment implemented in the simulation. This is important to judge upon 
the reliability of the conclusions derived from the latter. 

For our study we started from the well proven software package GEANT 
(version 3.21 distributed with CERNlib version 99) written at CERN [64]. In 
particular this frame comprises very powerful routines for the description of 
the detector geometry (only the CDC, Barrel and the magnet, return iron 
yoke, CDC/Barrel support structures and target passive volumes have been 
implemented for the ongoing attempt). Once the material of each volume 
and its characteristics have been specified the standard package routines are 
capable of calculating the particle transport in the magnetic field throughout 
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the detector including the decay with the appropriate hfetimes and branching 
ratios and the interactions with the penetrated media. Each decay product 
can then in turn be transported and let decay until they are all outside the 
interesting regions. If a volume is specified as sensitive (e.g. the CDC gas or 
the Barrel scintillator) hits are generated with the corresponding positions 
and energy-loss values. 

The remaining parts under the user responsibility (except the geometrical and 
material descriptions) are the event generator and the hit processing routine 
necessary to build the final output. 

With the thermal source a class of simulations, indicated in the following 
as "single-0", has been run. However the realistic track multiplicity of the 
full nuclear Ni+Ni event, originating the 0-meson, will certainly reduce the 
efficiency. In fact, due to the low granularity of the Barrel, one or more of 
the 180 scintillator bars can be crossed by multiple particles in the same 
event. In such sole improper position is recovered and both hits are 

discarded in the CDC/Barrel matching phase. If one of them corresponds 
to a kaon from the decay, the itself is lost. The CDC behavior too is 
not expected to be sampled correctly under the simplified situation with just 
two kaons not involving noise hits together with intersecting tracks and their 
mirror images. Each of the preceding reasons of inefficiency is affected by 
the distribution of the emitted 0-meson azimuthal angle with respect to the 
reaction plane where trajectories tend to concentrate. To take into account as 
much as possible the effect of the environment surrounding the paths of the K'^ 
and /^"-mesons created by the decay, an other class of simulations, named 
in the following as "embedded-0" , was accomplished. In it, to the thermally 
generated a whole IQMD Ni+Ni central collision (6 = 0) was superimposed. 
The orientation of the reaction plane, with respect to the (p azimuthal emission 
angle, was changed at random. We remind that the version of IQMD [44] 
employed contains no kaons, but only pions, nucleons and heavier fragments. 
In the following we will concentrate on single-0 Monte-Carlo where the simpler 
condition allows a detailed analysis of the stability of the code adopted, a clean 
comparison with the data and a reasonable estimate of the systematic errors. 
The embedded-^ simulation is cited to quantitatively evaluate the additional 
inefficiency coming from the nuclear environment surrounding the (p. 

In the hit processing routines a description of the front-end electronics has 
been embodied with the associated noise. For the CDC to derive from the in- 
formation provided by GEANT (position and energy-loss) the contents of the 
100 Mhz Flash ADC memories, four steps are necessary. First, for all tracks in 
a sector the charge corresponding to the given energy-loss is distributed with 
fluctuations along the path within each individual drift cell between electron- 
ion clusters. The latter are afterwards propagated under the position depen- 
dent influence of the electric and magnetic flelds [54] up to the sense wire. For 
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every cluster a generic pulse shape at each of the ends is also delivered. Second, 
all the pulses arriving at one specific readout channel are added together time 
bin wise without regarding the particle that produced them (in such way a 
cross-talk between different tracks in a sector is introduced). Third, the noise 
is added in an uncorrelated manner, separately for each wire, and in a corre- 
lated one, injecting charge simultaneously into all wires of one sector (as it was 
found to occur in the experiment). In both Gaussian distribution is 

assumed. Finally, the sum is digitized with a 8 bit nonlinear (10 bits effective) 
ADC accuracy. For the Barrel the time and energy digitized information is de- 
rived on the basis of a constant propagation velocity and attenuation length 
in the scintillator bars and light guides. A Gaussian smearing is added both in 
time and pulse height to include experimental resolutions. DoTible hits in the 
same scintillator bar are treated in a fashion that the fastest signal reaching 
the photomultiplier defines the time tag (independently for both ends) and 
the pulse heights are summed. Now the events are in a form compatible to 
be fed as input to the tracking and reconstruction code with which data are 
elaborated. 



3.2.1 The basic input resolutions as compared to the experiment 

Since the characteristics of the detection system can not be probed by a well 
defined beam hitting into the acceptance of the CDC, the performance of the 
Monte-Carlo digitization routines has to be monitored and tuned from the 
response to usual tracks. For those, of course, the real transverse momentum 
Pt or i?ia,b angle are not known. To remain as close as possible to the basic 
physical resolutions connected to how the chamber operates we have decided 
to compare the mean rms of the track fit residuals in the transverse plane 
RMSa;y, the mean rms of the track fit residuals in the 2;-coordinate RMS^ and 
the mean rms of the truncated energy-loss distribution RMSdE/dx- The results 
are displayed in fig. 6 for K'^ and X^-mesons from the single-0 simulation 
and K'^ alone from the Ni-|-Ni collisions due to hmitations in K~ statistics. 

For the data a set of identified kaons, defined as in the 0-meson analysis, was 
selected (see table 2). It has to be stressed that at a fixed p the energy- loss 
for unlike particles is not the same and the CDC behavior is correspondingly 
changed; so not all of them together, but exclusively ir"^-mesons, should be 
compared with the calculation. 

In the case of the Monte-Carlo only kaons known to have reached the Barrel 
sensitive volumes during the GEANT transport phase and with generated mo- 
mentum less than 0.6 GeV/c were considered. Because at most two tracks per 
event are present, the position of the primary collision vertex in the transverse 
xy-pla.ne can not be pinpointed as in the data; rather the refit is always forced 
to the origin where the initial position of all the reaction products is located. 



14 



As default do and zq are taken from the free fit. Such procedure, adopted here 
and in what follows, is expected to reproduce best the experimental situation 
where the number of CDC tracks is far greater 30) and hence the deter- 
mination of the event vertex xy-coordinates is accurate. Note that the true 
conditions correspond more to the embedded-0 scenario. 

To gain statistics a 130 MeV thermal source was used and the kaon decay 
inhibited. A comparison with a temperature of 70 MeV or with the kaon 
decay enabled has exhibited negligible differences below the systematic ones 
between data and simulation (this is the case also for the observables studied 
in the next section). 

The agreement is quite satisfactory for RMSa;^^ and RMSdE/dx (see fig. 6 a and 
c) . There are indications that the increase of RMS^ with pt is too weak in the 
simulation (see fig. 6 b). The stronger rise in the real detector is probably to be 
ascribed to a form of common high frequency noise induced in all wires, known 
to be relevant in this type of chambers [56], but not handled properly in the 
CDC digitization routine developed for GEANT. In fact while its repercussions 
are relatively negligible for RMSa;^, being the transverse position deduced from 
the drift time, they are rather important for RMS^ where the coordinate is 
determined from the charge division. 

The discrepancies between the values of RMS^^, and RMS^ in the Monte-Carlo 
and the data are small with respect to the much larger cuts A^yh^ and A^hu 
controlling the tracker behavior and so no readjustment in particular of the 
latter was done for the analysis of the simulation (see table 2). 

Once RMS^^y and RMS^ have been fixed, the CDC accuracy in transverse 
momentum RMS^^ and in t^iab angle RMS^ can be investigated with the 
Monte-Carlo incorporating the description of the detector response. Energy- 
loss, energy-loss straggling and multiple scattering in the passive and active 
volumes of the CDC/Barrel combination were accounted profiting of the de- 
scription of the physical processes embodied in the GEANT software pack- 
age [64]. The results are summarized in fig. 7 both for the mean systematic 
shifts and the resolution of the reconstructed quantities. 

For Pt a maximum systematic variation of ±1% is seen. The low pt distortion 
in fig. 7 a is to be attributed to the energy-loss in the target and chamber en- 
trance materials, as is proven by the simulation where these passive volumes 
were removed. The dissimilar behavior between positive and negative particles 
in the high pf region (see fig. 7 a) is due to the sectors tilting that introduces 
a systematic asymmetry. In fig. 7 b the pt resolution is also compared with 
what is expected from the analytical formulas of Gluckstern [65,66] consider- 
ing the multiple scattering in the CDC gas (dashed line) and the transverse 
hit position measurement error RMSj;^^ of 400 /im, as deduced from fig. 6 a. 
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and under the vertex constraint (dot-dashed hne). The qiiadratic addition of 
the two contributions (continuous hne) is quite near to what is obtained from 
the Monte-Carlo environment without target and chamber entrance materials, 
after tracking. It can be concluded that the multiple scattering in the chamber 
gas as well as in the passive volumes is the dominant limitation in the extrac- 
tion of transverse momentum through the entire pt range of interest. In the 
real experiment, if trajectories need to be merged between adjacent sectors, 
the situation is expected to be worsened by other problems, not included in 
the digitization routines. 

For the i^iab angle, the systematic effects can be neglected as compared to the 
corresponding accuracy. The worsening of RMS^ toward 90° can be understood 
as a 2;-coordinate constant resolution (not function of z itself) which translates 
into an uncertainty in -diah angle increasing toward ^?iab = 90°. 

The performance of the Barrel is less critical to model as it contains basically 
nothing else but the time resolution (being the position determined by the 
time difference) and the outcome will be supplied in the next section in terms 
of the final quantities. 



3.2.2 The shapes of the distributions associated to the selection cuts 

The comparison and tuning of the basic resolutions RMS^^y, RMS^ and RMSdE/d: 
is a prerequisite, but for the estimation of the detection efficiency the more 
complex derived observables, necessary for kaon selection in the 0-meson pro- 
duction evaluation, have to be studied: do, zq, CDC/Barrel matching in z- 
coordinate and in (/9-angle A^p, Cmass and Bmass. For each of these, 
attention should be devoted to both the mean and the width. In all cases the 
K'^ and K~ from the single-0 and the K'^ from the Ni-|-Ni FOPI data have 
been chosen. In the simulation, again the sample of particles known to have 
reached the Barrel sensitive volumes during the GEANT transport phase was 
always taken (even for variables like do) and the initial momentum was limited 
to 0.6 GeV/c. The values of the data refer to kaons identified with the cuts 
of the analysis (see table 2). Only the constraint relevant for the quantity 
under consideration was removed to see the shape of the distribution. A Gaus- 
sian fit was found satisfactory in all cases except the CDC-Barrel matching in 
</9-angle, where the influence of the scintillator bar finite width is present. Let 
us discuss all of them in turn. 

For do and Zq the results of the comparison are given in fig. 8 a and b, respec- 
tively. The mean and width of d^ are in good agreement, the sole discrepancy 
is in the bump starting at —0.5 cm in the simulation. The feature is connected 
to the problem that parts of the mirror tracks are sometimes also misinter- 
preted as real new ones. In the Monte-Carlo it happens in 0.07% and 3% of 



16 



the cases for the images of and i^'^-mcsons, respectively (with a statis- 
tics of 75000 tracks). The aheady mentioned tihing of the sectors renders the 
situation asymmetric and the missidentification is not equal for positive and 
negative particles. The bump is not visible in the data for a simple reason of 
statistics: the numbers of and K~ in the acceptance are the same when 
exclusively the (/)-meson decay is assumed as a source, however in the exper- 
iment the first is greater than the second by 2 orders of magnitude [39] and 
so the incomplete rejection of the image of negative kaon trajectories at the 
level of a few percent contributes a 10~^ background for positive ones. 

The width of the calculated zq distribution is about 4 times as wide as mea- 
sured. The reason is an underestimation of the correlated electronic noise 
pick-up that adds pulse height coherently into all preamplifier channels on 
one side of the CDC. Due to the projective geometry of the chamber this 
causes a strong distortion of the angle while the intercept with the beam axis 
is relatively stable. Since we operate with very wide windows on zq no attempt 
was made to fine tune its behavior. The cut at ±10 cm in the true conditions 
does not discard any good track but reduces the background from fake ones 
and from trajectories originated by secondary vertices, correspondingly the 
permitted range in the processing of GEANT output was opened to ±20 cm. 

The expected shift of the mean as well as the precision for the CDC-Barrel 
difference in both ^-coordinate and (/9-angle (Ai^ and A^, respectively) are 
gathered together with those of the data in fig. 9. For A^, the systematic 
displacement from of the latter (see fig. 9 a) is motivated by the difficulty of 
calibrating the longitudinal coordinate in the real experiment, however it is less 
important than the corresponding standard deviation by at least a factor of 
^ 2. The accuracy in A^ is somewhat better in the Monte-Carlo (see fig. 9 b), 
so the limits of the allowed region in its analysis were rescaled in the ratio 
6/8 from ±25 cm to ±19 cm (see table 2). In the simulation, the systematic 
oscillating shift in (y9-angle matching is imputable to a small inconsistency in 
the description of the CDC geometry, it is anyhow 3 times less than the 
resolution. The rms of /\^p is dominated by the shadow of the scintillator bar 
physical size and is hence easily reproduced. 

The comparisons of mean position and width of the Cmass and Bmass distri- 
butions are shown in fig. 10 in momentum bins of 0.1 CeV/c up to 0.6 CeV/c. 
As dimensionless fractions Cmass /itlk — 1 and Bmass /itlk — 1, where mx 
is the nominal charged kaon mass, were preferred to the bare Cmass and 
Bmass, respectively. With this choice both the mean and the standard devi- 
ation (7 can be expressed in percentage. All trends are similar. For Cmass the 
systematic shift (see fig. 10 a) reflects how well the employed empirical param- 
eterization is adapted to the energy-loss versus momentum dependence. It is 
a function of the calibration constants and was tailored to follow as closely as 
possible the experimental one (obtained in a far more complex environment 
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with many particle species and hence subject to more restrictions than in the 
Monte-Carlo). A residual uncertainty at the highest momenta is possible be- 
tween the 10% of the simulation and the 5% that could be extrapolated from 
the data with a straight line. The width in Cmass (see fig. 10 b) is dominated 
by the RMS^E/da;, found to be in good accordance (see fig. 6 c), it is then not 
surprising that the measurements are appropriately modeled. However sub- 
stantial extrapolation of the experimental points is necessary to compare the 
resolution to the simulation at the highest allowed momenta and as a conse- 
quence a margin of variation between the 8% constant value of the first and 
the 13% of the second is possible. For Bmass less freedom is available since 
no empirical parameterization is needed and the discrepancy in mean position 
(see fig. 10 c) could not be reconciliated, it is anyhow at worst of 2%, much 
less than the associated precision. Its accuracy is dominated by the scintillator 
Barrel time of fiight resolution (found to be 400 ps), again good agreement 
exists with the data (see fig. 10 d), but substantial extension is needed at the 
highest momentum and a discordance remains from the 10% predicted at 0.6 
GeV/c and the 8% that could be extrapolated with a straight line from the 
latter. 

The estimated possible maximum discrepancies in mean or width (for all rele- 
vant cases the distributions are Gaussian both in the data and the simulation 
so that the standard deviation a can be adopted) discussed in detail in this 
paragraph are summarized in table B.l. They will be the base to estimate the 
associated systematic uncertainty in the extraction of Cdet (see appendix B). 

The accepted ranges in the analysis of the GEANT and real results are col- 
lected in table 2. 

In the Monte-Carlo the contribution of each cut can be sorted out within the 
statistical accuracy of ±1%. The exclusion of hit multiplicities lower then 30 
causes no decrease and a reduction of 97% on the negative and positive track 
numbers, respectively, owing to the more frequent rejection of mirror im- 
ages, that has already been discussed in connection with the distribution 
shape. The efficiency of the do, selection itself, for the same problem, is also 
charge asymmetric being 95% for positive and 98% for negative tracks. The 
Zq and Cmass allowed windows do not leave out anything. The requirements 
on and are strongly correlated, so that once the (/9-angles selection is 
imposed, not much more is eliminated by the one on the 2;-coordinate difi^er- 
ence. They will be lumped together as a single condition in the following. The 
match, Bmass and momentum cut all have an efficiency near 97 — 98%. The 
previous three figures have then to be squared to take into account the posi- 
tive and negative kaon identification in every event to know the impact on the 
pair (i.e. 0) population. The 3% missing after matching can not be ascribed 
to the A<^ and A^ cuts, that correspond to more than 3 cr, but rather to the 
scintillator efficiency. Although for a path impinging along the normal to the 
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bar entrance surface, the depth of 4.05 cm is more than enough to generate 
a hght signal above threshold in 100% of the cases, it is nevertheless not true 
if the trajectories are curved by the magnetic field, and edge effects become 
important. It should be noticed that the width is 3.25 cm, to be compared 
with the already mentioned depth. 

A contamination at the level of (1.3 ± 0.2)% of the i^'~-couples with both 
momenta below 0.6 GeV/c from the part of the kaon spectrum above 0.6 
GeV/c is found with the present resolution of the CDC in the simulation. 

3.3 Final efficiency 

The position and width of the invariant mass peak are determined 

only by the pt and -i^iab angle systematic shifts and resolutions because, as 
already mentioned, once a particle is selected as a kaon candidate the nominal 
mass is assumed. If the peaks are fitted with Gaussian in the region of interest 
(1.01 — 1.03 GeV, as utihzed for extracting the number of 0-mesons in sec. 2.3) 
the mean is found to be 1.022 ± 0.003 GcV (data) and 1.020 ± 0.003 GcV 
(simulation) and the standard deviation a to be 5.8 ±1.4 McV (data) and 
4.7 ± 0.3 (simulation). However Lorentzian tails outside the fit region are 
distinguishable solely in the Monte-Carlo; the experimental spectrum is too 
much altered by the lack of statistics and the presence of the background. 
To grasp an idea of how looks the comparison to the 0-meson intrinsic V 
of 4.43 MeV, the FWHM of the Gaussians can be used: they are a factor 
of a; 3 larger. The agreement between the real and expected invariant mass 
resolution confirms that the broadening should be attributed to the CDC 
momentum uncertainty. Since the first and second moments are compatible in 
the measurements and the simulation, for calculating the fraction of identified 
(f) an identical binning is performed and the counts in the same 4 consecutive 
channels are added together. From the Monte-Carlo it can be estimated that 
(10 ± 3)% of the events are lost in such a way. 

Adopting the previous procedure the last two columns of table 3 have been 
obtained. The detection efficiency of the CDC/Barrel subsystem of FOPI for 
the 0, independently of the temperature in the considered range within one 
standard deviation, amounts to 

edet = 70% ± 1% . (7) 

The total efficiency with a thermal source of a temperature of 130 MeV is 
expected to be 

etot = (0.95 ± 0.02)% (8) 
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and with a temperature of 70 MeV to be 



etot = (0.24 ±0.01)%. 



(9) 



Up to now a single-0 has been used to derive the figures supphed in eqs. (7) 
through (9), but to really arrive at the final ones the additional bearing of the 
whole nuclear event surrounding the 0-meson has to be included. The lack in 
CDC/Barrel matching, due to the double hits in the latter, can be directly 
estimated from the data. From the simulation the expected Barrel detection 
efficiency for one charged kaon in the appropriate polar angle acceptance is 
(84 ± 2)%, owing to the azimuthal coverage and edge effects. It is difficult 
to compare safely with the experimental corresponding one, because and 
can be cleanly tagged only with Bmass and because of the low statistics. 
We opt to refer to a sample of well sorted protons in the CDC, in a similar 
momentum range of 0.4 — 0.6 GeV/c, although the energy-loss is not alike. The 
measured value is then of (85 ±2)% for few-track events, in good agreement 
with the prediction. On the other hand, for multiplicities, typically associated 
with 0-candidates, it drops to (81 ± 2)%. We introduce a correction to edet 
given by the square of their ratio 



This is just a lower limit on the reduction in the reconstructed fraction in- 
duced by the presence of all the collision products, it is nevertheless very well 
compatible to the embedded-0 GEANT outcome of (92 ±5)% obtained with 
a thermal source of at a temperature of 130 MeV superimposed on IQMD 
Ni+Ni events [61]. Such agreement suggests that double hits in the Barrel are 
the main origin of failure caused by a realistic multiplicity. Instead changes in 
behavior of the CDC tracking code seem to be minor, confirming its reliability. 

The maximum systematic error for the extraction of the efficiency from the 
simulation is estimated to be 40% comprehensive of all contributions: different 
behavior of each cut in the Monte-Carlo environment and in the data, small 
discrepancies in the geometry of the real and the simulated detector and the 
influence of the nuclear event surrounding the (see appendix B) . 



4 0-meson production yield 

We are now in a position to deduce the (/)-meson yield in central Ni+Ni colli- 
sions. Let us focus on the phase space portion relevant for the detector subsys- 
tem (see table 1) to initially avoid the uncertainty in the source temperature. 
Using the number of ^-candidates of 23 (see sec. 2.3) obtained from 4.7 • 10^ 



e^trk = (91±3)%. 



(10) 
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central events (see sec. 2.2), the 70% detection efficiency (see eq. (7)) di- 
minished under the reahstic track multiphcity (see eq. (10)) and taking into 
account the partial width of the into the iir~-channel of 49.1%, we find 
the 0-production probability in the CDC/Barrel acceptance per central event 
(0 < 6 <Ri 3.3 fm) Pace = (1-6 ± 0.5 ± 0.8) • 10"^ The first error is statisti- 
cal and the second systematic, calculated with the quadratic and linear sum, 
respectively, of the contributions of the data (see sec. 2.3) and the simulation 
(see sec. 3.3 and appendix B). 

An extrapolation to a full azimuthal angle coverage of the Barrel and no edge 
effects can be performed with e* (see table 1 and appendix A): this benefits 
that it is independent from the chosen source temperature 

P2. = Pacc/e* = (1-9 ± 0.6 ± 0.95) • 10"^ , (11) 

where the errors have an equal meaning as in the previous Pace- Paying at- 
tention to the kaon decay in fiight (as described in appendix A), P2tx can be 
immediately compared with theoretical predictions within the 'j?iab angle and 
total momentum acceptances. 

If the assumption of a source temperature of 130 MeV is made (see sec. 3.1), 
for the yield estimated in the full solid angle, it is found 

-^47r @ 130 MeV ~ -Pacc/emax®130MeV — (1-2 ± 0.4 ± 0.6) • 10~^ . (12) 

Finally for a source temperature of 70 MeV (see sec. 3.1), the same quantity 
of eq. (12) amounts instead to 

■^47r @ 70 MeV ^ -Pace/ emax@70MeV = (4.5 ± 1.4 ± 2.2) • 10~^ . (13) 

The errors in eqs. (12) and (13) have been derived like in eq. (11). 

The very preliminary result announced earlier in [49] is, when compared to 
eq. (12), smaller by a factor of ^ 2.5: this has to be ascribed to the previous 
version of the Monte-Carlo, less accurate than the current one and not yet 
tuned through a careful comparison with the real response as in sec. 3.2.1 
and 3.2.2. We underline that the actual size of candidates sample has not 
been changed with respect to both [49] and [61]. It is the extracted efficiency 
that was different, demonstrating the importance of a sound study about the 
agreement between simulation and experiment for the observables necessary 
to sort out the kaons and showing how appropriate it is to estimate a max- 
imum possible systematic errors (see appendix B) as attempted in eq. (11) 
through (13). 
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The (t>/K~ ratio extrapolated to the full solid angle can now be given. The 
/^"-mesons production probability in 47r can be deduced to be (2.1 ±0.4) -10^^, 
using the rapidity distribution measured by KaoS [38] also for ^^Ni+^^Ni at 
1.93 AGeV and for central collisions, but with the KaoS experiment centrahty 
selection (0 < 6 < 4.4 fm) and assuming a Gaussian shape. Scaling this value 
with the mean number of participants as Apfi-t"'^ [67], the corresponding one to 
the FOPI experiment central trigger is (2.6 ± 0.5) ■ 10^^. If the last procedure 
is applied to /T^-mesons consistency is found between the KaoS data [38] 
and those of FOPI [37]. The 0//s:--ratio is then (0.44 ± 0.16 ± 0.22) and 
(1.7 ± 0.6 ± 0.85) at 130 MeV and 70 MeV temperature, respectively. Hence 
referring to the highest, but still reasonable of the two, and minding at the 
branching ratio, it can be established that at least ~ 20% of the /('"-mesons 
arise from the decay of an intermediate 0. 



5 Compcirison with theoretical models 

The first detailed theoretical study of 0- meson yield in the Ni+Ni system at 
1.93 GeV was that of Chung, Li and Ko [68,69,70]. Their relativistic transport 
model RVUU describes the rapidity distributions of p and tt data and the 
collective flow reasonably well. For the 0-meson production, the calculation 
foresees four scenarios: i) no in-medium effects, ii) in-medium effects only on 
the kaon masses, iii) in-medium effects on the kaon masses and through this on 
the decay width and iv) in-medium effects also on the 0-meson mass. The pa- 
rameterization of [28] is taken for the in-medium modification of the 0-meson 
mass together with the elementary cross sections from their own model [68]. 
Within this frame the dominant channel for ^-production is txN (j)N . The 
sensitivity was assessed relative to the free scenario i) for the i^~-meson to 
be a factor of 3 more in case ii), iii) and iv) and for the to be a factor of 
2 reduction in case iii) and 10% less in case iv). Their expectation for the 
0-meson occurrence probability in central collisions is 2.3 ■ 10"^ in case iv). It 
is extracted from [69] with the appropriate mean over h (eq. (5)). Assuming a 
(pN elastic cross section of 0.56 mb (their preferred choice) and of 8.3 mb (the 
maximum compatible with the (/)-meson photo-production data) the deduced 
final kinetic energy distribution exhibits an inverse slope of 110 and 130 
MeV, respectively. Adopting a temperature of 110 MeV for extrapolation and 
comparison to the experiment, the model underpredicts the data by 6 times. 
Such discrepancy is, however, not too surprising when viewing the large un- 
certainties in the parameterization, close to threshold, of the elementary cross 
sections. 

Recently Barz, Zetenyi, Wolf and Kampfer [71,72] have been engaged in a 
new effort, profiting of the measurements now available for the elementary 
pp — > (f)pp reaction at ^/s — y^Sth = 83 MeV [48]. It is based mainly on 
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their previous work [73] about the role of elementary three body interactions, 
where the intermediate particles can be off-shell. It is found that in current 
microscopic transport models a small error is done considering successive two 
body collisions, where every particle is on-shell, if all the relevant intermediate 
steps are included. They point at the Bp and 7rp-channels as important ones, 
neglected in any preceding investigation, and develop their own predictions 
starting from this remark. The in-medium effects are accounted for the 
alone with the same result of [28] as in the calculations exposed before [69]. 
The (t)N elastic cross section is placed at 0.5 mb, near to the preferred estimate 
of [69] . The awaited 0-production probabihty for central collisions (9% of the 
total cross section) is 1.7 • 10~'^. The authors apply the FOPI CDC/Barrel 
geometrical acceptance filter and obtain a value of 2.7 • 10^^ in reasonable 
agreement with the experimental finding of 1.9 • 10^^ (see eq. (11))- It is 
however not clear how well the implementation elaborated for the moment, 
can explain the global experimental observables. In addition it remains to be 
seen how in-medium modifications of the kaon properties change the picture, 
since the yield is coupled to the sum of the and K~ masses through its 
decay width. More work is needed to compare carefully their theory with the 
experimental achievement reported here. 



A completely different approach is that of the thermal models [74,75,76] where, 
assuming a chemical equilibrium before freeze-out, the abundances can be 
predicted (not the yields themselves). Particular care is to be exercised in 
the treatment of conservation laws like baryon number, electric charge and 
strangeness. In the thermodynamic limit both canonical and grand canonical 
formulations are equivalent, but for small systems, however, the discrepancies 
are large. In fact in general all these quantum numbers should be handled 
in the canonical treatment (implying exact conservation) as has been proven 
necessary for strangeness, while for baryon number and electric charge the 
grand canonical one (implying conservation in mean) is adequate [74]. If also 
strangeness is described with the latter, the data in the SIS energy regime are 
not matched. Determining the temperature and baryonic chemical potential 
from the d/p, t:^ /p and j'n" relative experimental intensities, the maxi- 
mum of the (^jK" one, that can be accommodated, is of the order of 0.1 [76]. 
The last value is compatible with the figures of sec. 4, within the large sys- 
tematic and statistical errors, only if a temperature of at least 130 MeV is 
hypothesized. This situation illustrates the limitations of the actual small ac- 
ceptance measurement where the source is not well characterized and of this 
class of models where the species kinetic temperatures are poorly predictable. 
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6 Conclusions 



A first experimental production probabifity of sub-threshold 0-meson in heavy 
ion collisions has been obtained. Before solely and K~ data were available. 
With a pseudo-vector meson a new sector of the theoretical machinery, pro- 
posed to describe in-medium effects, has been tested. It supports the emerging 
awareness of a strong interconnection between the , K~ and 0-channels. 

Here only the data gained from the CDC/Barrel subsystem are discussed; 
the effort necessary to make a quantitative statement for the probability of 
occurrence in central collisions is described. The creation of a 0-meson is 
a very rare event at sub-threshold energies: out of 4.7 • 10^ events in the 
sample of central collisions, barely 23 candidates were found after background 
subtraction. Despite the poor statistics, a signal to background ratio of ~ 1 
could be achieved when applying appropriate cuts. 

An important part of the effort has been the estimation of the efficiency that 
was splitted into two parts: one, emaxi due to the geometrical acceptance, al- 
lowed momentum range and and decay in flight and a second, edet, 
originating in the identification requirements necessary to decrease the back- 
ground. While Cmax is between 1.35% and 0.36% (depending on the source 
temperature), showing that the greatest loss in events is to be attributed 
to the localization of the phase space covered by the detector around target 
rapidities, e^ct is 70% indicating the validity of the cuts chosen in reducing the 
background without rejecting too many good events. 

Furthermore edet was proven to be insensitive, within the statistical error of 
the simulations, to the source temperature (and it could be argued even from 
the shape of a moderately anisotropic one); on the contrary e^ax varies by a 
factor of 4 when the temperature changes from 70 to 130 MeV. The comparison 
with calculations or new results is hence greatly simplified: in fact while emax 
can easily be evaluated along the guidelines of appendix A by other groups, 
edet can exclusively be extracted from a complete simulation of the detector 
response, requiring a specific knowledge of the apparatus; it has been deduced 
once for all in the present work. 

Future confrontations will benefit from a thorough discussion of the possi- 
ble sources of systematic error: a value of 40% and 9% has been given for 
simulation and data, respectively. 

This measurement suggests that, near threshold, a large fraction of the K~- 
mesons (at least 20%) arise from the decay of an intermediate (/)-meson. As a 
consequence it would be important to verify, both experimentally and theo- 
retically, to what extent the K~ phase space distribution and the 0/iir~-ratio 
can be contaminated from the 0-production channel. The discrepancies with 
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the predictions based on microscopic transport codes can partially be ascribed 
to insufficient knowledge of the elementary cross sections. A possible role of 
intermediate p-mesons production for the 4> yield appears to be an interesting 
aspect of investigation. The (f)/K~-Ta,tio seems to be hardly compatible with 
thermal models. 

In order to improve the knowledge on the experimental side, a more complete 
coverage of the acceptance is needed. Especially, the access to the region of 
low transverse momenta in the proximity of mid-rapidity is important for a re- 
liable reproduction of the total yield of such rare particles produced in central 
heavy-ion collisions. A first step along this way is taken by the ongoing effort 
to include the Helitr on/Forward Wall data [50]. In a near future the FOPI 
collaborations plans to accumulate more data with an increased statistics on 
0-production, when the undertaken upgrade program will be finished. On the 
other hand, the HADES collaboration will be able to detect the 0-meson 
through the leptonic channel in search for direct evidences of in-medium ef- 
fects. Both plans make the subject open for promising new developments. 
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Table 1 

Geometrical acceptance of the CDC and the Barrel. 
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Table 2 

Momentum, tracker and identification cuts used for the analysis of the data and of 
the simulation. 
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eBar (%) 
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Table 3 

Geometrical ecDC and eear, maximum emax, total etot and detection e^et (see eq. (6)) efficiencies for four thermal source temperatures. 
The geometrical contribution is given for the CDC alone (ecDc) and the CDC/Barrel combination (eear)- The latter €Bar includes the 
not complete azimuthal coverage of the Barrel and edge effects (see e* of table 1 and appendix A). If, within the acceptance of the 
CDC/Barrcl subsystem (eBar), the upper momentum limit at 0.6 GcV/c is imposed, Cmax is obtained. For eBar and Cmax the reduction 
due to the kaon decay in flight during their path to the Barrel is considered. All figures assume that the only open channel for the 
0-meson is the K~ one. 
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Fig. 1. Upper panels: logarithm of the truncated mean energy-loss dE/dx versus 
momentum p/Z for positively (panel a) and for negatively (panel b) charged parti- 
cles reaching the Barrel. The dashed lines represent the empirical parameterizations 
used for the energy- loss versus momentum in the Cmass reconstruction algorithm; 
the constants chosen are the same as those of the data analysis for the extraction 
of the (/>-candidates. Lower panels: momentum p/Z versus velocity v for positively 
(panel c) and for negatively (panel d) charged particles. The dashed lines represent 
the relativistic momentum versus velocity relation with the nominal masses. 
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Fig. 2. Cmass versus Bmass for positively (panel a) and negatively (panel b) 
charged particles with momentum less then 0.6 GeV/c. The area enclosed by the 
dashed line corresponds to the kaon identification cuts used in the <j)-meson analysis 
(see table 2). 
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Fig. 3. Charged particle multiplicity (PMUL) distributions corresponding to the 
Minimum Bias trigger (empty) and the Central trigger (shaded). The upper scale 
in impact parameter b is derived with the assumption of a sharp cut-ofF geometrical 
model. 
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Fig. 4. Invariant mass spectra for identified pairs with the final cuts (see 

table 2). The histogram is the signal; the two lines enclose one standard deviation 
from the background reconstructed with the event-mixing technique when the mix- 
ing stack is filled with Ni-|-Ni collisions containing either a or a Er~-meson; 
the error bars refer to the background reconstructed when the mixing stack is filled 
with reactions containing both a and a K~ . 



32 



-2 -1.5 -1 -0.5 0.5 1 

y 




Fig. 5. Transverse momentum over mass pt/m versus rapidity y plots for the orig- 
inal </>-meson (panel b) and for the and i('~-mesons produced in its decay 
(panel a). The GEANT simulation with a source temperature of 130 MeV and in- 
cluding the decay in flight of the kaons is used. The rapidity is renormalized such 
that —1 is the target, -|-1 is the projectile rapidity and is mid-rapidity, where the 
(/> source is located. The continuous and dashed lines in panel a and b illustrate 
the acceptance limits of the CDC/Barrel combination (compare with table 1). The 
dot-dashed curves in panel b represent the most probable momenta for a relativis- 
tic Maxwell-Boltzmann distribution of emitted (p centered at y = and with a 
temperature of 130, 110, 90 and 70 MeV, respectively. 
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Fig. 6. Results obtained from the CDC digitization routine for the mean rms RMSa;^ 
of the track fit residuals in the transverse plane (panel a), the mean rms RMS^ of 
the track fit residuals in the z-coordinate (panel b) and the mean rms YIMS^^;/^^ 
of the truncated energy- loss distribution (panel c). Note that in panel b and 
if~-mesons from the simulation lie on top of each other. The range in energy-loss 
of the data (500 — 1600 a.u.) corresponds to kaons with momenta 0.1 < p < 0.6 
GeV/c. 
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Fig. 7. Transverse momentum and longitudinal angle T?iab systematic shift and res- 
olution from the Monte-Carlo. In addition to K'^ and i^~-mesons from the standard 
environment used in the present work, also from a simulation without target 
and chamber entrance materials are considered for comparison. In panel b the "pi 
root mean square deviation is compared with what is expected from Gluckstern's 
formulas [65] describing the effect of multiple scattering in the CDC gas (dashed 
line) and of a transverse hit position uncertainty of 400 /xm (see fig. 6 a) under the 
vertex constraint (dot-dashed line). The two contributions are added quadratically 
(continuous line). 
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Fig. 8. Comparison of distributions in do (panel a) and zq (panel b) of i^+-meson 
tracks. The solid histograms are from the data and the dashed ones from the simu- 
lation. In panel b a Gaussian fit is also shown. 
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Fig. 9. Comparison of the means and widths for CDC/Barrel matching in 
z-coordinate and ip-angle (A^ and A^, respectively) between and ivr~-mesons 
in the simulation and in the data. 
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Fig. 10. Comparison of means and widths of the non dimensional quantities 
Cmass/niK — 1 and Bmass/rriK — 1 between and i^~-mesons in the simu- 
lation and in the data. 
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A Direct calculation of the geometrical acceptance 



If the 0-meson source is not spherical or the temperature is outside the range of 
table 3 (so that an interpolation can can not be used) Cmax has to be calculated 
from the events generated by the assumed model. We briefly describe the filter 
procedure, with the related setup parameters, which is needed to compare 
future data or the prediction of a microscopic transport code with the 0- 
production probability measured here in the CDC/Barrel acceptance. After 
the (f) has decayed into a K'^ and a X^-meson when both have a momentum 
smaller than 0.6 GeV/c, for each of them it is necessary to determine if it 
hits the nominal Barrel surface. The latter is defined as a cylinder located at 
a radius i?Bar with a z extension from z^i,^ to z^^^. First for the inverse of 
the radius of curvature in the transverse plane A;, that can be obtained from 
eq. (1), the condition 2/k > R^ax has to be fulfilled. In this case the length It 
of the transverse projection of the track from the origin to the Barrel can be 
calculated as 

and the 2;-coordinate of the hit in the scintillator bar is 

z^hk, (A.2) 



where h is the helix pitch defined hy h — Pz/pt- Second the value given by (A.2) 
has to lie inside the interval Zmin < z < Zmax with the z-adn ^-nd ^^max of table 1. 
To account for the possible decay in fiight, each kaon must be weighted with 
a survival probability 

e.p L^^] , (A.3) 



where ttik, tk, Pk are the kaon mass, mean life and momentum, respectively 
and I is the corresponding distance of fiight to the Barrel 



/ - /* Vl + ¥ . (A.4) 



After the weighting with the product of the quantities (A.3) for each particle 
of the couple, the final outcome delivers , where 



"-max 



(A.5) 
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In order to obtain emax of eq. (6) and table 3, a reduction factor e*, compre- 
hensive of both kaons, has to be apphcd to compensate for the incomplete 
azimuthal coverage of the Barrel due to the regions occupied by the CDC 
support structures and to the inter-module and inter-bar spacings. Because 
the depth of the scintillator bar is bigger than the width, edge effects are 
important at the level of a few percent (sec sec. 3.2.2) and have to be in- 
cluded as well. The value of e* has been estimated through the ratio between 
the output of the filter just outlined and of GEANT where the geometry is 
defined accurately implementing separately each scintillator bar with its ap- 
propriate positions and dimensions. It turns out from the simulation that e* 
is independent from the source temperature in the considered range, within 
one standard deviation, and can be determined once for all as (81 ± 2.5)% 
(see table 1). In a comparison of theory to the experiment one can benefit of 
eq. (11) and barely supply so that 

P4. = P2./e^3x ■ (A.6) 



B Estimate of the systematic errors 

Our estimate of the efficiency is biased by systematic errors of distinct types: 
a non-identical behavior of every cut in the data and simulation, small geo- 
metrical discrepancies between the real detector and the simulated one, the 
infiuence of the nuclear event surrounding the 0-meson. We discuss each of 
them {do, zq, Bmass, Cmass, A^, A^, nutsi upper momentum, Barrel geom- 
etry, multi-track environment) in turn. 

Beginning with the K~-m.eson selection requirements, the ones on do and 
zo are assumed to contribute negligible systematic distortions. As a matter of 
fact, the distribution of the first is well reproduced and the bump starting at 
—0.5 cm in the simulation is due to mirror tracks that have to be discarded 
(see sec. 3.2.2). Any eventual correction arising from the improper purge of 
such cases is assumed minor. The latter does not reject good events and is 
wide enough to accommodate the small discordances in the experimental and 
awaited shapes, once the condition imposed on zo in the analysis of GEANT 
is released as given in table 2. 

The Bmass, Cmass and observables obey a Normal law with good ap- 
proximation in the data and in the simulation, besides such assumption is 
legitimate for the estimation of systematic errors. As long as the real and 
calculated spectra have a Gaussian shape there could be only two sources of 
concern about the effect of a fixed accepted window: a difference in the mean 
or in the standard deviation a. In truth if both are not quite similar the inte- 
gral inside the same range is not equal, as can be easily obtained from the error 
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function. The discrepancy in the efficiency of a cut between the simulation and 
the data translates directly by twice the amount (once for each kaon of the 
pair) to the correctly identified fraction e^et- The change for variations in the 
mean or the a judged reasonable from our previous discussion (see sec. 3.2.2) 
is written up in the rightmost column of table B.l in percentage of the number 
of events allowed by the selection. Since the domain of uncertainty quoted in 
mean or a are the maximum expected ones, an analogous statement is true 
also for the systematic error (in the same spirit for evaluating the effect of 
the shift in mean always the larger of the two cited cr's is preferred). In the 
case of Cmass the permitted region is widely asymmetric and solely the lower 
exclusion value matters, this peculiarity was accounted. 

For A^, the width of the required interval used in the analysis of the data 
was rescaled for the simulation one. A small residual repercussion of the sys- 
tematic discrepancy in the standard deviation is expected as compared to the 
others where the window is left untouched, here we will neglect it. Even if 
A(p, partially reflecting the geometrical shadow of the scintillator bar, does 
not properly follow a Normal law, a Gaussian is not a too bad approximation 
for estimating the systematic error. It turns out to be very small, after all 
the cut is sufficiently wide, and it can be neglected (see table B.l). The two 
rms of A^ are well compatible with each other and no adjoined systematic 
distortions are considered from this side. The alteration in rejection caused 
by the dissimilarities in the shape of real and foreseen A^ and Ac/, distribu- 
tions is in general not expected to be the major contribution to the systematic 
error in the matching ratio. In fact the corresponding accepted windows are 
~ 3(7 wide and, as already discussed (see sec. 3.2.2), the origin of the loss is 
rather the scintillator efficiency coupled to the edge effects for curved tracks. 
However, if the matching fraction for single kaons predicted by the GEANT 
model (84 ± 2)% and the measured one for well defined CDC protons tracks 
(85±2)% are compared (see sec. 3.3), the discordance of 1% is not statistically 
meaningful. So no further contribution is introduced from the systematic lim- 
itations in the implementation of all the detailed processes affecting A^ and 
A^. 

The invariant mass spectrum, being dominated by the CDC momentum res- 
olution, is also approximately Gaussian (see sec. 3.3). No statistically signif- 
icant disagreement exists between measurements and Monte-Carlo either for 
the mean or for the width, so the associated systematic error is assumed to 
be negligible. 

The probability for a track to have an assigned number of hits does of course 
not follow a Normal law but is rather asymmetric with a tail towards low hit 
multiplicities (even for kaons reaching the Barrel) populated by fake tracks 
alone in the simulation and also by pion background in the experiment. It is 
difficult to quantify the reliability of the description of this quantity because it 
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depends on many subtle particularities of the detector response, so an estimate 
(see table B.l) equal to the remainder of what is rejected in the simulation 
(3%) and in the data (5.5%) was taken, it can not be exceed by the true 
discrepancy. 

We interpret systematic errors as reasonable estimates of the possible displace- 
ment of the true result by effects not completely controllable and opt for a 
linear sum with the aim of arriving at trustable gap in which the final effi- 
ciency should lie. If the disagreement with a calculation is outside this, serious 
doubts should be shed on the theory itself. Adding all the discussed sources 
a maximum systematic error of 13% is found for one kaon and of 26% for the 
pair, i.e. for edet- 

The margin of uncertainty due to the upper momentum limitation was esti- 
mated varying the cut of ±1%, that is roughly the systematic distortion of 
the momentum reconstruction (see fig. 7) and taking the biggest change (see 
table B.l) in the size of the selected 0-sample (the source temperature was 130 
MeV). Including this last contribution gives an expectation of 28% systematic 
error imputable to imperfect characterization of the kaon selection constraints 
behavior. 

The difficulty in the determination of the geometrical acceptance can be re- 
duced to the problem of knowing with enough precision the respective posi- 
tions of the target and the Barrel, whose polar angle coverage in the forward 
direction is smaller than that of the CDC. Since the target is replaced after a 
set of runs, its displacement has to be deduced each time from relative mea- 
surements of the parts of its holder system, a maximum safety boundary of ±1 
cm is appropriate. A simulation was run with a thermal source of 0-mesons 
at 130 MeV temperature and the Barrel was moved of ±1 cm, the influence 
on Cmax is reported in table B.l. 

For estimating the possible systematic mismatch in the reduction of detection 
efficiency by the multi-track environment, the difference in the total percentage 
of correctly identified events in the single-(/) and embedded-0 Monte-Carlo is 
indicated in the rightmost raw of table B.l. 

A maximum comprehensive systematic error for the extraction of the efficiency 
from the simulation of 40% is estimated. 
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Table B.l 

Contributions to the systematic error. The values of the selection windows considered are those of the analysis of the data and the 
simulation (see table 2). For the first four cuts the mean and standard deviation a quoted for the data are based on estimations of 
the maximum discrepancy and for the simulation are simply what is obtained at the highest permitted momentum of 0.6 GeV/c. From 

the difference between the measurements and the Monte-Carlo the systematic error on the cut effect is calculated assuming a Normal 
distribution and given as a percentage of the events allowed. For the last four contributions the estimation is adapted for each as explained 
in the text (see appendix B). 
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